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Abstract 

The vortex dynamics of Euler's equations for a constant density fluid 
flow in is studied. Most of the paper focuses on singular Dirac delta 
distributions of the vorticity two-form uj in M^. These distributions are sup- 
ported on two-dimensional surfaces termed membranes and are the analogs 
of vortex filaments in and point vortices in M^. The self-induced velocity 
field of a membrane is shown to be unbounded and is regularized using a 
local induction approximation (LIA). The regularized self-induced velocity 
field is then shown to be proportional to the mean curvature vector field of 
the membrane but rotated by 90° in the plane of normals. Next, the Hamilto- 
nian membrane model is presented. The symplectic structure for this model 
is derived from a general formula for vorticity distributions due to Marsden 
and Weinstein [20|. Finally, the dynamics of the four-form w A a; is exam- 
ined. It is shown that Ertel's vorticity theorem in M^, for the constant density 
case, can be viewed as a special case of the dynamics of this four-form. 
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1 Introduction. 

The Navier-Stokes and Euler equations for a classical fluid in a physical domain of 
spatial dimensions greater than three have typically attracted mathematical inves- 
tigations aimed at understanding fundamental issues related to the equations and 
their solutions. On an n-dimensional Riemannian manifold, for example, these 
equations have been investigated for their deep geometric, topological and func- 
tional analytic properties — good introductions to such research with many impor- 
tant references can be found in the books by Arnold and Khesin [6J and Temam 
ll32l . There is also a slender body of theoretical and computational work — the 
latter no doubt facilitated by major advances in computational power in the last 
decade — on the Navier-Stokes equations in higher dimensions with the focus on 
turbulence. A few examples of such work are [l3Tlf23l[TTl[T3l[T0|. 

It is fair to say that, in general, classical fluid flows governed by these equa- 
tions in higher dimensions-even if one restricts to M'^-have been far less investi- 
gated than flows in M.^ and M^. In particular, little attention seems to have been 
paid to the development of low-(phase-space)-dimensional mathematical models 
of flows in M^. Motivated by this, this paper focuses on a set of mathematical mod- 
els of vorticity, in the framework of Euler's equations, which have found much 
favor amongst mathematicians, physicists and engineers for studying flows in 
and R^. These are models of the dynamics and evolution of coherent vortices in 
which the vorticity field is a singular Dirac delta field. Examples of such models 
are point vortices, vortex filaments/rings and vortex sheets. A good introduction 
to these models can be found in Saffman [27|. 

The importance and significance of vorticity for flows in three or less dimen- 
sions is so well-established that it needs no elaboration. It is therefore to be 
expected that it plays an equally important role for flows in higher dimensions. 
Vortices, as is well-known, also occur in other areas of physics most notably in 
superfluids. For this paper, however, the focus will be entirely on classical fluids. 

The outline of the paper is as follows. In section 2, some basic facts of the vor- 
ticity two-form in are presented. In section 3, the notion of a vortex membrane 
in is introduced. A membrane is a singular Dirac delta distribution of vorticity 
supported on a two-dimensional surface. It is shown that the self-induced velocity 
field of a membrane is unbounded and a local induction approximation (LIA) is 
used to regularize the velocity. In section 4, the Hamiltonian membrane model is 
introduced. It is shown that the evolution equation for a membrane — again, af- 
ter a LIA regularization — has Hamiltonian structure relative to a symplectic form 
derived from a general formula for vorticity distributions due to Marsden and We- 
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instein [|20ll . In section 5, the relation of the (regularized) velocity field to the mean 
curvature vector field of the membrane is made clear. In section 6, the dynamics 
of the four-form oj A u — which is identically zero in and lower dimensions — 
is examined and it is shown that Ertel's vorticity theorem in M?, for the constant 
density case, can be viewed as a special case of the evolution equation of this 
four-form. 



2 Basics of vorticity in E^. 

In M^, the vorticity field is a genuinely two-form field. It can no longer be as- 
sociated with a vector field, as in M^, or with a scalar field, as in M?, but can 
be represented, if necessary, as a skew- symmetric matrix field. It follows that a 
natural setting for investigating vorticity in higher dimensions is the exterior al- 
gebraic/differential geometric framework, along with the associated Hamiltonian 
formalism, pioneered by Arnold [l4l|5]| and developed, among others, by Marsden 
and Weinstein [i20il . 

Equipping with the Euclidean metric, the 1-form associated with the ve- 
locity field (through the metric) v{xi,t) = {vi{xi,t)), i = 1, ■ ■ ■ , 4, is 

4 

= ^^fjfixj (1) 

i=l 

satisfying the divergence-free condition 

I)iv{v) = -dv^ = *d*v^ = 0, (2) 

where 5 : f]'=+i(M^) n''{R'^) is the co-differential operator and * : fi^(R^) 
^4-fc^]^4^ is the Hodge star operator^ Here, ^l'' represents the space of differential 
/c-forms on M^. 

The vorticity two-form is 

u = dv^ (3) 

/ dvi dv2\ , , / dvi dv3\ , , / dvi dv4\ , , 
= -7^ \- t: — dxi A dx2 + -t: h -7: — dxi A dx^ + h 7^ — dxi A 0x4 

\ 0X2 0X1 J \ 0X3 0X1 J \ 0X4 OX\ J 

( dv2 dv3\ , , / dv2 dv4\ , , / dv3 dv4\ , , 

+ -7^ \- 7^ — dx2 A dx3 + -7- h 7^ — dx2 A dx4 + h 7^ — dx3 A 0x4 

\ 0x3 0x2 J \ 0x4 0x2 J \ 0x4 0x3 J 

=: uji2dxi A dx2 + cuisC^Xi A dx3 + w^fixi A dx4 

+ UJ23dX2 A dx3 + UJ24dX2 A dX4 + +UJ34dX3 A dX4 (4) 
'For the minus sign in front of 9, see 121, page 458 
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Closedness of w. Since uj is exact, it is closed i.e dw = 0. However, there seems 
to be no obvious way of identifying du = with a divergence-free condition as 
in M"^ or M^. In the above coordinates, the condition gives rise to the equation 

^^^^^ dx, Adx2Adxs+(^-^ + ^) dx, A dx2 A dx. 







dxs 


0X2 




duO] 


\ (9X4 


dx, 



duji2 


duji4 


dx4 


8X2 


f duj23 




\ (9X4 


dx. 



OX\ ) \ 0X4 0x2 0X1 

(5) 

+ ( „ + -TT—] dxi A dxs A dx4 + ( — ^ + | dx2 A (ixa A dx4 

0x1 J \ 0x4 0x3 0x2 J 

(6) 

at all points in M^. This gives rise to four equations which can be written in 
compact form as 

1^ = 0, (^, J, k) = {(1, 2, 3), (1, 2, 4), (1, 3, 4), (2, 3, 4)} (7) 

with the convention that ojki = —coik if k > i and so on. 
The Lie-Poisson equation for vorticity evolution [20] is 

-^ + C,LJ = (8) 

Using exterior algebraic properties of Lie derivatives [El, it can be shown that the 
equations of each component of cj — in coordinates — are 

^+ E ^+ E ^ = (9) 

In the above form of writing the equations, the convention adopted is ujkj = —^jk 
if k > j and Uik = —oom if i > k. 



Vortex surfaces. The essentially two-form nature of cu in MJ^ means that there 
are no vortex lines, i.e. curves tangent to vorticity vector fields [|27ll since the latter 
do not exist. From an exterior algebraic viewpoint, a vorticity vector field in is 
the vector field {*ujY, where * is the Hodge star operator and (j is the operator that 
relates 1-forms to vector fields via the metric. In M"^, however, *uj is a 2-form and 
there is no canonical way of identifying {*uj)'^ with a vector field. 

To elaborate, consider, more generally, (local) submanifolds 5 C M" to which 
the two-form u is perpendicular or, equivalently, 

i*suj = 0, (10) 

where : 5 — t- M" is the inclusion map and superscript * denotes puUback of 
forms. Since * * a; = (— l)^*^"^^^a;, the above equation is equivalent to saying that 
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*uj is tangent to S (see BH, p.540, for an introduction to perpendicular and tangent 
forms). For *u to be a non-trivial form on S, it follows that the dimension of S 
must be equal to or greater than the order of the form i.e. dim(iS) > n — 2. The 
lowest possible dimension of S is therefore one@in and two in W^. Hence, in 
there are no vortex lines — instead, one thinks of vortex surfaces which satisfy (flOl) 
at each point. 

3 Vortex membranes and 3-branes. 

Singular Dirac delta vorticity fields: vortex membranes and 3-branes Re- 
call that in M?, when the voriticy is a singular Dirac delta field supported on a 
S with (Dim(5) = 0, Dim(5) = 1) one obtains the (point vortex, vortex sheet) 
models. Similarly, in M'^, when the voriticy is a singular Dirac delta field sup- 
ported on a iS with (Dim(iS) = 1, Dim(iS) = 2) one obtains the (vortex filament, 
vortex sheet) modelsH 

Singular Dirac delta vorticity fields in that are the analogs of point vortices 
and vortex filaments in and MP, respectively, are two-dimensional surfaces 
S — henceforth referred to as membranes — such that 

i^u = 0, 

at all points of S, where z : S — )■ is the inclusion. The precise sense in which 
the vorticity field is a Dirac delta field supported on S is made clear further down. 

Note that a vortex sheet in M^, when embedded in M^, is not a vortex mem- 
brane because the vorticity two-form acts on different planes for a sheet and a 
membrane. For a sheet embedded in W^, the vorticity two-form at any point acts 
in a plane that intersects the sheet non-trans versally in M^. But for a membrane 
the vorticity two-form at any point acts in a plane that intersects the membrane 
transversally in M^. 

Just as vortex filaments in can be 'stacked together' in such a manner that 
the vorticity two-form continues to be perpendicular to the resultant surface, to 
form a vortex sheet, similarly, membranes in can be stacked together so that 
the vorticity two-form continues to be perpendicular to the resultant volume to 
give analogs of vortex sheets which are termed 3-branes. 

For the rest of this paper, however, the focus will be only on membranes. 

^And zero in consistent with the fact that vorticity in the plane is equivalent to a scalar 
function. The above geometric argument also shows why 'vortons' in M!^ are not genuine repre- 
sentations of vorticity; for more on vortons see, for example, Novikov | |24) or Winckelmans and 
Leonard 1 36 1 

■'The models are nothing but special cases of the M'^ models and it is sometimes useful to 
adopt this viewpoint. 
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Vorticity two-form for a membrane. To write the vorticity two-form for a 
membrane in an appropriate basis, recall an elementary fact about 2D surfaces 
in M^. Unlike in where the normal vector is defined uniquely at any point of 
a 2D surface, there is no unique vector at each point of the surface but rather a 
plane of normal vectors. At any point p G S, the local decomposition thus exists 

TpR^ = TpS © iVp, 

where Np is the span of the normal vectors to S at p. And so the most natural way 
to write the vorticity two-form, which is perpendicular to S, is as a volume form 
on Np — similar to how Marsden and Weinstein EOl wrote it for vortex filaments 
in M'^. As they noted, the form ij/i is a volume form on the plane of normals at 
each point of a curve in W^, where fi is the volume form on M'^, t is the unit tangent 
vector to the curve and i is the interior product [l2l. To see this, simply choose an 
orthonormal basis (t, ni, for at each point of the curve, where (rii, 77-2) is a 
basis for the plane of normals, to obtain i^yU = dni A dn2 

Similarly, for membranes choose an orthonormal basis {ni,n2, ^1,^2) for 
at each p G S, where ^1,^2 G Np and (^1,^2) £ TpT.. The volume form on in 
this basi^ is 

/i = drii A dn2 A dti A dt2 (11) 

and one obtains 

ifji^j/i = drii A dn2 (12) 

as a volume form on Np. 

It follows that the vorticity two-form for a single vortex membrane of strength 
r in can be written as 

uj(m) = TS{m — (i^j/i), m G M}, p G S, 

= T6{m — p)dni A dn2, (13) 

where 5{m — p) is the Dirac delta function supported on S. 

Note that membranes in W^, like filaments in M'^, are boundaryless due to 
the closedness of the vorticty two-form. Co-dimension two vortex surfaces, such 
as membranes, have also been examined in the context of the Ginzburg-Landau 
equations by, for example, Lin [fT9l and Jian and Liu [fTTll . 

''The volume form can also be written in the canonical basis /i = dxi A dx2 A dx^ A dx4 and 
noting that dxi A dx2 A ^0:3 A dx4 — dni A dn2 A dti A dt2 since the two bases are related by an 
orthogonal transformation. 
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Assumption. For the rest of the paper, the membranes considered will be as- 
sumed to be compact and arcwise-connected subsets of W^. The former assump- 
tion is made simply to avoid technicalities of Stokes' theorem for non-compact 
manifolds and it is conjectured that most of the results that follow also hold for 
non-compact membranes. The latter assumption is needed for the proof in Ap- 
pendix |Al From a physical point of view, neither of these assumptions are consid- 
ered restrictive. 

Strength of a membrane. Before proceeding further, the notion of the strength 
of a membrane, F, needs to be defined and clarified. 

Definition 3.1 The strength of a vortex membrane S is defined as 



where A is a smooth two-dimensional disc-like surface which intersects S transver- 
sally at a single point, such that dA (a closed curve) cannot be shrunk to a point 
without intersecting S, and ia '■ A is the inclusion map. Equivalently, using 

Stokes theorem and since u = dv\ 



1. Since the intersection is transversal, S, which is of co-dimension two, inter- 
sects A, also of co-dimension two, only at isolated points. 

2. As in R'^, the definition using the notion of circulation, (fT4l) . can be taken 
as the basic definition of F and the other definition, using vorticity, can be 
viewed as a consequence of Stokes theorem. 

3. Since is a Dirac delta field, the integral in the above definition should be 
interpreted — like in — as the limiting value of the integral defined for a 
tubular neighborhood of the membrane. This is made clear below. 

4. For a given OA there can be many As but, as in M^, Stokes theorem applied 
to A or, equivalently, (fT4l) shows that the choice of A is irrelevant. 

The next step is to show that the above definition is independent of the choice of 
dA. The proof of this is presented in Appendix lAl 





(14) 



where iqa '■ dA — >■ A is the inclusion and iqa := ° idA- 



Notes: 
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3.1 Self-induced velocity field of a membrane. 



In this subsection, the self-induced velocity field of a membrane is examined. It is 
shown that — as for filaments in — the field diverges and requires a regulariza- 
tion. A local induction approximation (LIA) is then used to show that the leading 
order term in the regularized field is proportional to the mean curvature vector of 
S rotated by 90 degrees in the plane of normals at each point. 

To show this, first express the velocity field of a given smooth vorticity field. 
The relationship V^^; = — V x a; generalizes in as 

So; = 9dv'', 

where 9 is the co-differential operator. The above is obtained simply by applying 
the co-differential operator to both sides of the definition of the vorticity two- 
form. The operator 9d on functions is nothing but the Laplacian in M^. Writing 

= ^11=1 '^i ^ ^^i' where : M"^ ^ M, obtain 

4 

9a; = 9d^^(t'j A dxi), 
1=1 

= {V'vf. 

The last line can be obtained by either expanding the terms on the right side of 
the previous line in coordinates or using properties of the exterior algebra objects 
(wedge product, d, 9 and *). Note that use is made of the divergence-free property 
of V. The equation is now equivalent to the Poisson equation 

V^Vj = *{dxi A *9a;), i = 1, • • • ,4 

Denoting by G{m, mo) the Green's function of the Laplacian in W^, where m, mo e 
R^, the solution to the above is 

Vi{m) = / G(m, mo) A *{dxi A *9a;) /x^q, i = 1, • • • , 4 

= -/ {{ei,{*{dG{m,mo) A*uj)f))i^^„ i = 4 (15) 

using properties of the exterior algebra objects and integration by parts. To apply 
the latter, it is of course assumed that u decays sufficiently rapidly at infinity 
for the boundary term in integration by parts to vanish. In the above, is the 
canonical volume form on where the subscript mo means that the variable of 
integration is mg, and ((;,;)) is the standard metric on R^. Note that the variable 
in the d operation is again mo. 
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To obtain the self-induced velocity field of S, substitute (|T3I) in (|T5l) and switch 
the basis from the canonical to (ni, n2, ti, ^2) to get 

vsi{m) = -T 6{m-p){*{dG{m,p) A{dtiAdt2))fl2p, meM.'^,peJ: 

S{m - p){—dn2 - —dnifnp, 



= r / 5(m - p)((VG ■ n2)ni - (VG ■ ni)n2)^p, 
= r J {{VG ■ n2)ni - {VG ■ ni)n2)up, m,peT. 
In the above, is the same as (fTTI) and z/p = dti A (it2- 



(16) 



Remark. There is no canonical way of fixing ni and 722 or ti and t2- The ve- 
locity expression should however be invariant to the choice of the normals and 
the tangents. But this is a straightforward check. The integrand in (fT6l) and 
Up are invariant under the (orthogonal) transformation in the plane of normals 
(ni, 722, ti, ^2) — > '"'2, ^1, ^2) and the (orthogonal) transformation in the tan- 
gent plane (ni,ra2,ti,t2) {ni,n2,ii,i2). 



Recall that in M^, 



G{m,p) = — —-, (17) 

r — r 1^ 



where r and r are the position vectors of points m and p, respectively. And so, 

7r"^(| r — r pj 



(gradient with respect to f). Substituting (1181) in (1161) . it is obvious that the in- 
tegrand develops a singularity as p — )■ m. The integral (lT6l) therefore needs 
to be evaluated as an improper integral. More specifically, the behavior of the 
leading order term of the integral as p — m is needed. Since there is no diver- 
gence in the integral due to non-local effects (as can be easily seen by noting that 
(VG ■ ni)njiyp = 0(l/f) as f — )■ 00), it suffices to evaluate the integral in a 
neighborhood of m. 

Consider two closed neighborhoods Us, f/e C S of m, such that C Us- The 
parameters 5 and e measure the size of these sets (discussed more below). Let 
U := Us\U^ denote the resulting annular domain and 

J := r / ((VG ■ n2)ni - (VG ■ n^)n2)up 
Ju 
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Proceed now to examine the limit 



Leading order term of ^^/(m) := lim I, 



(19) 



with Us fixed. 

Coordinates. All the integrals presented above are integrals of differential forms 
on smooth manifolds and by their definition such integrals are evaluated using co- 
ordinate charts or, in other words, mapping the domain of integration smoothly 
to/from R". The value of the integral is of course independent of the choice of the 
coordinates but a judicious choice simplifies the evaluation. In particular, a set of 
coordinates with the following two properties are desirable: (i) the level curves of 
the coordinates intersect orthogonally and (ii) the coordinates denote arc-lengths. 

Consider the following system of curvilinear polar coordinates consisting of 
smooth curves parametrized by arclength si measured (in the metric) from m 
and all passing through m. Call these 'diameter' curves and these do not meet 
except at m. Each of them connects two points on the boundary of Us- Similarly, 
consider a set of closed curves enclosing m parametrized by arclength S2 such 
that each closed curve intersects every diameter curve orthogonally at exactly the 
same length of the diameter curve from m. To each point in Us is thus associated 
a unique pair of coordinates (si, S2) £ and this defines a coordinate map : 
(si, S2) p. Define also a coordinate 9 — the angle made between the tangent 
vector to a 'radial' curve at m and any fixed line in T^S — such that S2 = /(si, 0), 
where / is a smooth (well-defined) function with /(si, 0) = 0. ^ is associated 
with the entire 'diameter' curve. 

And so / is evaluated as 



The puUback of the integrand is simply the integrand composed with 0. To eval- 
uate 4)*{i'p), first note that since the choice of (ti, ^2) is arbitrary, fix these vectors 
by taking them to be tangent to the level curves of si and 9 at each point of Us- 
Next, write down the coordinate map explicitly using the coordinates of p 



and a straightforward computation, using properties (i) and (ii) of the coordinates, 
gives = dsi A ds2 = dsids2- And so 

/■/(5,27r) p5 



= r / / ((VG • 712)711 - (VG • 7ii)7i2)(0(si, /(si, e)))g{s^, 9)ds,d9, 




: (si, S2) ^ {x{si, S2),y{si, S2), z{si, S2), o(si, S2)) 



(20) 




(21) 
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where g{si,9) = df{si,9)/d9. Since the level curves of si shrink to zero as 
Si — )■ 0, one observes that g{si,6) — )■ as si — )■ (for all 9). 

At each point in Us, there exists an orthogonal frame 9),n2{si, 9),ti{si,9),t2{si,9)) 

and, in particular, at the point m, there exists the one-parameter set of frames 
(?2i(0),n2(0),ti(0,6'),t2(0,6')) tangent to the 'diameter' curves. Wlog, letti(0, 0) = 
ti and ti{0, n/2) = ^2- Taking r = 0, expand f for small Si and fixed 9 along each 
'diameter' curve 

Similarly, for i = 1,2, 
raj(si,6') 

- «.(o) + ..||(o,.) . 1^(0,9) + 1^(0.9) + 1^(0.9) + 0(,5), 



and 



For small si, /(si, 9) ~ Si6' since the mesh of coordinate curves must approach a 
mesh of polar coordinates in the tangent plane at m. And so set ^5^(0, 9)/dsi = 1 
in the above series expansion. 

Using (fTSi) . the integrand of / (modulo the constant factor 2T/n'^) is 

(f • n2)ni — if ■ ni)n2 



f 1^ 



Substituting the series expansions for small si, (suppressing the (51,6*) notation 
temporarily) 

1 1 



1 A . , , ^/ 2^ 



Next, 



(r ■ n2)ni — (r ■ ni)n2 



dn2 1 , , dti \ , . ( drii 1 , , dti , , , 
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so that (re-introducing the {si,6) notation), 
(f ■ n2)ni — (f ■ ni)n2 



1 



f 1^ 



rr. n^ ^^2(0) 1 (9tl(O,0)\ 



ri2(0) 



0(1), 



1 



iMo,o).^)-.(o) 



iMo,«).^l,.(o) 



+ 0(1), 



using the fact that ti(s, ^^) ■ r2i(si, 
ating the limit (fT9l) one obtains 

Leading order term of vsi{m) 



V(si, 6). Substituting in and evalu- 



lim 



>0 TT^ 



loge 



27r r 



9^2(0) 

9s 1 



rii(O)- ti(0, 



9ni(0) 
9s 1 



ri2(0) 



(22) 



In conclusion, the self-induced velocity field of a membrane has a logarithmic 
divergence, as do vortex filaments in [|27l |3l. The direction of this infinite 
velocity is determined by the integral in the above expression and this is evaluated 
in a later section. 



Regularization of the infinite self-induced velocity using the LIA. As for fil- 
aments in M^, the infinite self- induced velocity can be regularized by a local in- 
duction approximation (LIA) [|27l l3l. This means that in the above evaluation one 
ignores the neighborhood of m in evaluating the integral. Equivalently, fix e at 
a finite, non-zero, small value. The resulting leading order term in vsi{rn) is then 
taken as the regularized self-induced velocity. And so 



vsi. 



reg 



m] 



log 



2iT r 



dn2{0) 
dsi 



t,(0,^).^)MO) 



dsi 



(23) 



e fixed, e << 6 



Note the slightly different way in which the above is written, compared to (1221) . 
The 0(1) term log S is retained to avoid attaching any spurious importance to the 
absolute value of e. Moreover, the use of the LIA and other regularizations for 
filaments in have shown that, from a physical point of view, the argument for 
the logarithmic term is best expressed as the ratio of two length scales, e and 5, 
with 5 representing a typical length scale of the vortical structure [27] . 
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4 Hamiltonian structure for the membrane model. 



In this section the Hamiltonian structure of the membrane model is presented. 
The Hamiltonian structure follows from the general symplectic/Poisson structure 
governing the dynamics of a vorticity field (in the framework of Euler's equa- 
tions). Marsden and Weinstein Il20l . expanding on the pioneering work of Arnold, 
showed formally that the appropriate phase space for the dynamics of a vorticity 
field specified by a vorticity two form is an infinite-dimensional manifold known 
as the coadjoint orbit of the dual of the Lie algebra g* (identified with the space 
of vorticity two-forms). The Lie algebra g (identified with the space of diver- 
gence free velocity fields) is associated with the Lie group of volume preserving 
diffeomorphisms of the flow domain (in general, an n-dimensional Riemannian 
manifold M). This manifold is equipped with a symplectic structure, known in 
general as the Kirillov-Kostant-Soreau symplectic structure on coadjoint orbits, 
and is given by 



In the above, the vorticity two-form u is identified with a point on the coadjoint 
orbit and (— —C^uj) — again, two-forms — are identified with tangent vectors 
at that point. Here m, f G g are identified with any two divergence-free vector 
fields on M and £ denotes the Lie derivative. More details can be found in EOl fTl 



As applications of (|24)) . Marsden and Weinstein further showed that (|24)) re- 
covers the symplectic structure of the well-known point vortex model in and 
provides a symplectic structure for the vortex filament model in M'^. 

Symplectic phase space for membranes. Consider now the phase space V of 
membranes i.e. the images of (unparametrized) maps, modulo reparametrizations 
with the same image, S : — )■ M"^. An element of P corresponding to such a 
map will be denoted by p. 

To apply (|24l) to membranes, first evaluate the tangent vector to coadjoint or- 
bits term for a membrane. Start with a tubular neighborhood C of the geomet- 
ric surface S as in Appendix |Al It is useful to think of C as a normal bundle 
over S. Choose coordinates (/i, ^2, si, 52) where (si,S2) are coordinates for E 
and {hih) are coordinates in the plane of normals. Choose an orthogonal frame 
(ni, 712, ti, ^2) at each point of C, tangent everywhere to the coordinate curves, 
and such that (rii, n2) lie in the fiber. Without loss of generality, assume the vor- 
ticity two-form in C to be a; = udrii A dn2, where a) is a constant. In other 
words, the vorticity two-form is the same on any plane of normals. As C shrinks 
to S, a) —7- 00 (such that the strength F of C remains the same) and one ob- 
tains the membrane vorticity two-form. Write the velocity one-form on C as 




(24) 



|2D. 
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= vi{p)dni + V2{p)dn2 + V3{p)dsi + Vi{p)ds2 where p = (/i, ^2, si, S2) . Use 
Cartan's 'magic' formula [|2T|| . the closedness of cu and v divergence-free, to ob- 
tain 

CjjOJ = diyoo = (vi^ + -^2^ I dni A dn2 = (v ■ VnO;) dni A dn2, 
V oh dh J 

where V„ is the gradient operator in the plane of normals. Note that (by the con- 
struction of C) the term in parentheses is a Dirac delta function on the boundary 
of C. In the membrane limit, 

Cyio = {\vn \ S{m — p)) drii A dn2, m E M.'^, p E H 

where Vn = ^/vf + f|, to give the identification of C^oo with a normal vector field 
onS. 



Remark. The conclusion is unaltered if a more general representation of u is 
chosen since in the membrane limit only the drii A dn2 component of to, with a 
coefficient that is constant (though infinite!), survives. 



Substituting tu for membranes, (fT3] ). in the right hand side of (|24l) gives T driiA 
dn2 (u, v) u. But since the form drii A dn2 acts on only the normal components at 
any point of S, the symplectic form on V can be finally written as 

n (X, Y) (p) = r / dni A dn2 (m„, w„) v (25) 

In the above, X,Y E TpV and m„, t>„ are the corresponding normal vector fields 
onS. 

Hamiltonian function. The Hamiltonian function on V is the kinetic energy of 
the fluid flow due to the membrane but regularized using the LIA. 
To write the fluid flow kinetic energy, 

K.E. = - [ A *V\ (26) 

as a functional of membranes it is useful to write it first as a functional of the vor- 
ticity two-form. This requires the introduction of the equivalent of the divergence- 
free 'vector potential' vector field A in M!^ which, it may be recalled, is related to 
the velocity field v inM.^ hy v = V x A, and to the vorticity vector field uj in 
by Poisson's equation: 

V^A = -uj (27) 
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In M^, since vorticity is a two-fomi field, the vector potential is also a two-form 
field satisfying similar equations as above. More precisely, the following is true 
in R^. For V divergence-free, there exists a two-form A, satisfying dA = 0, such 
that 

-*d*A = dA = V^ (28) 

In the above, 9 is the co-differential and acting on 2-forms in (see [2J), p. 457). 
Apply the exterior derivative operator to both sides to obtain the equivalent of (|27l) . 

ddA = 00. 

The operator d5 is the Laplace-de Rham operator without the 9d part. The equa- 
tion is then equivalent to a Poisson equation for each component of the two-forms 
as shown below. Let 

4 4 

A = Aijdxi A dxj, 

i=l j>i 

where Aij : M"^ — )• R. Computing d9^ using = 0, obtain 



4 4 

dg^ = ~ X] X] "^"^^ijdxi A dXj, 

1=1 j>i 



which is equivalent to the following set of equations in Cartesian coordinates 

V^Ai, = -coij, J > % (29) 

Existence and uniqueness of A, for given w, therefore follows from the existence 
and uniqueness of solutions to Poisson's equation in M"^. The solution to (l29l) . 
using (fTTI) . is 

where, as before, r and f are position vectors of points m and p, respectively, and 
the subscript on fi denotes integration with respect to p. 

For a membrane, switching to the frame (^1,71,2,^1,^2) with origin at some 
p E H and using the notation 



A{m) = At-^t^{m)dti A dt2 + An^t-^{m)dn2 A dti + ■ ■ ■ , m G 



d4 



(with the correspondence 1 ^ ni, 2 ^ n2,3 ^ ti and 4 o ^2), equation (|29l) 
becomes 

V^A^,^,{m) = -T6{m-p), (30) 
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and the right side is zero for all other components of A. Note that the choice of p 
makes no difference since the right side is always a Dirac delta function supported 
on S. Replacing the notation A^^^j by the more meaningful notation As, obtain 

A^{m) = T [ —-^-—-dh A dt2, (31) 

with the other components of A taken, without loss of generality, to be equal to 
zero. 

Returning to the kinetic energy integral (|26|) . proceed to write it for mem- 
branes. Using the above results and integration by parts, first rewrite it as 

K.E. = 1 I V A *gA 



u A*A 



Substitute ^13 to obtain 



r2 



K.E = - [ As{m)u, m e E 

-Up u, m,p E (32) 




2 J^h7r^i\r-f\') 



Of course the integral is not convergent till one regularizes the self-induced con- 
tribution. Once appropriately regularized the above integral can be viewed as a 
functional of S — since it depends only on the strength of the membrane and its 
geometry. And this leads to the Hamiltonian model for membranes. 



Regularization of the kinetic energy using the LIA. The regularization of 
the kinetic energy using the local induction approximation (LIA) proceeds in a 
similar way to the regularization of the self-induced velocity using the LIA. In 
other words, one removes a neighborhood of m in evaluating the inner integral 
in (|32l) . This removes the blow-up in the inner integral and the outer integral is 
then evaluated as a regular integral. As in the regularization of the self-induced 
velocity, there is no blow-up due to non-local effects. To obtain the leading order 
term of the resulting kinetic energy integral for e small, it suffices to evaluate the 
inner integral on a subdomain of S-in particular, the annular domain introduced 
previously U := Us\U^. 

f f 1 

K.E.req '■= Leading order term in — / / —-- ttv^'p t « 5 

2 J^Ju \\ r -r\^) 
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Using the expansions previously computed for small si, with r = 0, obtain for the 
inner integral 



2tv rS 



dsid9, 



^'Jo J. (|f(«i,e)P) 

Jo Je [Sl 

-lege, + 0(1), 
TT 



In the above, 



S 

Ce := - >> 1 
e 



And so the regularized kinetic energy which is also the Hamiltonian function H : 

if(S) := K.E.reg ■= — log c, / (33) 

Hamiltonian vector field. To compute the Hamiltonian vector field correspond- 
ing to (l33l) and (l25l) . the first task is to compute the variational derivative of the 
Hamiltonian. 

Let 5r : S — !■ correspond to a variation (5S G TpV. The former is a field of 
normal vectors on S. Elements of T*V are also identified with normal vectors on 
S and are paired with elements of TpV using the Euclidean metric in i.e. 

(us, as) := y << Xu, X^ » z/, G TpP, G Tp*P, 

<< Xu,Xa » dsids2 

where Xu, Xa are the corresponding normal vector fields (note that this pairing 
is degenerate if not restricted to normal fields). The second integral implicitly 
assumes the existence of a partition of unity and a set of coordinate charts covering 
S. Here, f/s denotes the entire domain of integration in the coordinate charts. 
The variational derivative — a normal vector field corresponding to an element of 
T*V — is then defined in the usual way: 

^) = /" « 6r, Xsh/sj: » dsids2 := lim - [i7(S + e5S) - if(S)] 

(34) 
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Consider the local chart defined previously (|20l) 

: (Si,S2) {x{si, S2),yisi, S2), z{si, S2), o{si, S2)) = r{si,S2) 

where, (si, S2) are again arc-length coordinates corresponding to an orthogonal 
mesh (but not necessarily the same polar mesh introduced earlier despite the use 
of the same symbol S2). The puUback of v := dti A dt2 G A^(M^) under this chart 
is defined as 

<i)*{dtiAdt2){u,v) := rfti Adt2(D0(M),D0(t;)), u, G T<^-l(p)M^ 
= {uiV2 - ViU2){ti ■ ti)(t2 ■ 12) 
^ (f)*{dti A dt2) = (ti ■ ti)(t2 ■ t2)dsi A ds2, 

where mi, U2 are the Si, S2-components of u etc. and dsi A ds2 G A^(M^). The 
terms ti ■ ti and t2 ■ ^2 are retained in this form for now since, as shown below, 
they contribute to the variations . 

Now consider the following local chart for the perturbed position of the mem- 
brane 

:(si,S2) 

(x + e6x, y + e6y, z + e6z, o + e6o){si, S2) 
= (r + e5r){si,S2) 

A slightly longer puUback computation then shows that 

(f)*{dti A dt2){u,v) 

= iUiV2 - U2Vi) ((ti ■ ti)(t2 ■ t2) + eiiSh ■ h){t2 " ^2) + (^^2 " t2)(tl " tl)) + 0(6^)) 
= {{h ■ ti){t2 ■ t2) + e((5ti ■ ti){t2 ■ t2) + {5t2 ■ t2)ih ■ h)) + 0(e2)) dSi A dS2, 

The functional derivative SH/ST., is identified with a normal vector field on S. 
From (|34l) . obtain 



(((5r(0(si, S2)), X5j^/52(0(si, S2)))) dsids2 



— log Ce - lim 



dti Adt2- / (iti A dt2 

S+e5S is 



■ ti){t2 ■ t2) + {5t2 ■ t2){tl ■ ti)) dSidS2 

{5ti ■ ti + 5t2 ■ 12) dsids2, 

To manipulate the above integral, consider the following 1-form in on S, 

(5r ■ ti)dt2 



(35) 
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Since S is compact and boundaryless, from Stokes theorem 

[ d{{5r-t,)dt2) = (36) 



Now, in a local chart, 



OSi 

d {{6r ■ ti)dt2) = I ^d^'^'^d^i'^"') ^2))^^i^^2, (37) 



Using ti = dr/dsi,t2 = dr/ds2 and making the usual assumption of variations 
commuting with differentiation, i.e. 

dSr d6r 

oh = TT—, dt2 = -TT— 
OSi OS2 

it follows from ^ and ^ that 

dti 

{6ti ■ ti) dsids2 = — I (t^ — ■ 5r) dsids2, 



/ iSt2 ■ ^2) dsids2 = - [ ■ 5r) dsids2 



Substituting back in (|35] ). finally obtain 

Since, Xsh/st, has non-trivial components only in the normal directions, we can 
write 

where ni{p) , n2{p) denote the basis vectors in the plane of normals at p. Using 
ti(p) ■ rijip) = Vp, the components on the right (in the ni(p) — n2{p) plane) can 
also be written as 

'drii drii \ , . 



9s 1 dS2 



etc. Letting, 



1 / dux drii dn2 dn2 , , , 
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it follows that 



XsH/5T.{p) = log C,^{p) 



TT 



The Hamiltonian vector field Xh is obtained from the general form of Hamil- 
ton's equation on symplectic manifolds [EH: 

Using (|25l) and the above derived result for the variational derivative, this equation 
becomes 

r f dni A dn2 (Xh, 5r) u = f \\^^-> logQ^)) dsids2, 

is i[/ECK2 \ \ TT / / 

/ {{Sr,{XH ■ ni)n2 - {Xh ■ n2)ni)) dsids2 



5r, log ) ) dsids2, 

Since 5r has non-trivial components only in normal directions, this gives 

2r, 

Xh = R-7T/2 log c^n, 

TT 

where -R_vr/2 = ( ? ) ' Hamilton's equations for the membrane evolu- 



tion is 



aE 2T 

= Xh = ■ — log c,^ (39) 

at TT 



5 The mean curvature vector field. 

In this short section, it will be shown that in (l38l) is the mean curvature vector 
field of S. Moreover, it will be shown that the right hand side of (|39| ) is the 
same as the regularized self-induced velocity field (|23l) obtained from a direct 
computation. 

Recall that for a 2D surface S in M^, the second fundamental form at p G S is 
the quadratic form 

&{V){p) := - {{Bn{p) ■ V), V)) , Ve T^S, 

where {n,ti,t2) are the unit vectors of a moving orthogonal frame and the unit 
normal vector is viewed as a map n : S — )■ 5^ so that T)n{p) : TpS(~ M^) — )■ 
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Tn{p)S'^{^ M^) is the derivative map [!8l[33l. The mean curvature n is then defined 
as 

For a 2D surface S in M'', there is a second fundamental form associated with each 
of the normal directions, Ui, of the moving frame (rii, ^2, ^1,^2)- 

6i(V)(p) := - ((Dni(p) ■ V), V)) , ©2(V^)(p) := - ((Dn2(p) ■ V), V^)) , 

V e TpE, 

The mean curvature vector K is then defined as (see, for example, [9|) 

2 1 

^(p) = E 2 + ^p^'^'^p^ ^^^^ 

and one checks that (l38l) is the same as the above. 

To show that the right hand side of (|39| ) is the same as the regularized self- 
induced velocity field (|23l ) is equivalent to showing that 

1. r%i(0,e) . ^^d9 = -J^(m) ■ ni(0) (41) 
27r Jo dsi 

^ [ \i{0,e) ■ ^^dd = -^{m) ■ ^2(0) (42) 
Zn Jq dsi 

(reverting to the notation used in the self-induced velocity field section). How- 
ever, this is a straightforward exercise and for the case of can actually be 
found as an exercise problem in, for example, Thorpe's textbook fSB] (Exercises 
12.16 and 12.17). Begin by noting that the integrands on the left are nothing 
but -©i(ti(0,^))(m), i = 1,2. Since ti = ti(0,0) and ^2 = ti(0,7r/2), write 
ti{0, 9) = cosOti + sin 6*^2- Substituting for ti{0, 6) and performing the integra- 
tions, one obtains that the integrals are equal to (— l/2)(©j(ti) + &i{t2)){'m). 
Using (go]), proves (gT]) and (l42l) . 

It is interesting to contrast the dynamics of membranes, as given by (|39| ). with 
the dynamics of the vortex surfaces considered in [lT9l [TTl where the surfaces 
move parallel to the mean curvature vector. 

6 The dynamics oiuo /\uj and an application to Er- 
tel's vorticity theorem in M^. 

As an application of vortex dynamics in M'', it will be shown in this section that 
Ertel's vorticity theorem in M"^, for the constant density case, is simply a special 
case of the dynamics of the four-form u /\uj. h few general ideas related to this 
four-form are discussed first. 
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Basics of a; A w. All four-forms are identically zero in and lower dimensions. 
The four-form a; A in M"^ is given by 



cj^ = 2 — UJ13UJ24: + U14U23) dxi A dx2 A dx^ A dx^ (43) 

It is an exact form by definition 

oj Au = dv^ Au = d[v^ Au) (44) 
and by Stokes' theorem it follows that 

u Alo = i* {v^ Alo) (45) 

DcR4 JdD 

implying J^^ a; A a; = if a; has compact support or decays sufficiently rapidly at 
infinity. 

It is also straightforward to see from ([8]), using derivation properties of the 
operators involved, that 

^<i^ + £„(.A..) = (46) 
or 

Further, writing 

u Au = f2fi = f2A fi, (47) 

where /i is the canonical volume form on R!^, it follows using derivation properties 
and the fact that v is divergence-free, i.e C^fi = 0, that 

^ + CJ2 = 0, (48) 

A consequence of the general results of Khesin and Chekanov IfTSl for an n- 
dimensional Riemannian manifold is that the Euler's equations in admit an 
infinite number of integrals of the form J^^ /(/2)/^- Arnold and Khesin |l6l call /2 
the 'vorticity function' and discuss several of its properties, see also the discussion 
in ll35l . The dynamics of w A is also studied in the numerical and theoretical 
work of (13). Interestingly, the earliest work known to the author that studies 
equation (|48l) . and generalizations, is a relatively lesser-known work of Hill [fTSl . 
also cited by Truesdell [[34ll (p. 173). Indeed, Hill examined Euler's equations on 
an n-dimensional Euclidean space. 

Application to Ertel's vorticity theorem. First, consider the case when the 
flow is with respect to a non-rotating, inertial frame. Consider /(x, y, z, t) : — )■ 
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M. Let whe a divergence-free velocity field in satisfying Euler's equation for 
a constant-density fluid flow 

dw ^ 

— — h w ■ Vw = — V-, 

at p 

where the V operators are in M?. The following vector field in R"^ 

V ■■= (wj), 

is therefore also divergence-free. For v to satisfy Euler's equation it follows that 
/ must satisfy 

f + u..V/ = 

Euler's equation, written for the velocity 1-form 

dil' 

— + = dg, (7 : ^ M (49) 

ot 

where the d operator is in M^, is then equivalent to the equations 

- + ..v. = -v(^j+v, = v^ — 

^ + ^ . V/ = ^ 
dt do 

Euler's equation for w and the conservation law for / are both recovered by choos- 
ing g{x, y, z, o) = (w^ + p)/2 - p/p. 

The six components of the vorticity two-form are then 

UJ = (wi,d/), 

where tUi is the vorticity two-form dw^. Denoting by X^- the vector field in 
corresponding to Ui, i.e. X^^ = (*]R3Wj)* = (^23, -^is, ^12), obtain 

a; A a; = 2{X^^ ■ \/f)dx A dy Adz Ado 

And so from (|46l) and (|48] ). one obtains Ertel's Theorem for the potential vorticity 
Xuii ■ V/ in a non-rotating frame ||25l . 

^ (X^, ■ V/) = 0, 

where D/Dt = d/dt + wV. 
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To demonstrate the result when the flow is with respect to a constantly rotat- 
ing frame — the usual setting of Ertel's theorem — the Lie-Poisson formulation of 
Euler's equation in a rotating frame in M" needs to be briefly discussed first. Re- 
call [6l, l205 that in non-rotating, inertial frames in M", Euler's equations can be 
written as Hamilton's equations associated with a Lie-Poisson bracket on g* — the 
dual of the Lie algebra of the Lie group of volume preserving diffeomorphisms of 
— as follows 



d[v^] 5F 



dt ' 6[v^ 



SF SH 



(50) 



Here, v E g is a divergence-free velocity field in the inertial frame and [f ^] = 
dv^ G g* is the equivalence class of 1 -forms generated by f ~ f + dh, : — )■ 
M, and identified with the vorticity two-form of v. The right hand side of the above 
equation is the Lie-Poisson bracket on g*. The pairing ( , ) between elements of 
g and g* is given by ([t>''], := J^„ v^(w) /i and is non-degenerate due to the 
L^-orthogonality of gradient fields and divergence-free fields. The Hamiltonian 
function H : q* R is the flow kinetic energy H{[v'"]) = (1/2) J^^v^ A 
and is a well-defined function, i.e. independent of the choice of the representative 
element of [v'"] used in the integral, because of the L^-orthogonality. Computing 
the functional derivative gives 5H/S[v'"] = v and using Lie bracket [ ] properties 
one obtains Euler's equations as 

m 

from which one obtains both (|49l ) and ([8]). 

To apply this formalism to a constantly rotating frame in R", let f G g be 
a divergence-free velocity field of the flow observed in the rotating frame. For- 
mally, g can be viewed as the Lie algebra of the Lie group of volume preserving 
diffeomorphisms of M" that go to elements of SO{n) at infinity. More specifi- 
cally, each volume preserving diffeomorphism in this group has associated with it 
a unique element of SO(n). Equivalently, each Vr has associated with it a unique 
element of so{n) — the Lie algebra of SO{n). The uniqueness follows simply by 
implementing the infinity condition in the inertial frame which is that the volume 
preserving diffeomorphism goes to the identity map or, equivalently, the velocity 
field goes to zero. Note that the uniqueness of course does not extend the other 
way. 

The dual of the Lie algebra, g*, is identified with elements [v^] = dvl = Ur, 
where tUj- is the vorticity two-form of the flow observed in the rotating frame. The 
equivalence class is again generated by gradient fields in the rotating frame that 
go to zero at infinity. The pairing is then defined as 
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where Q^, Qu G 5o{n) are constant two-forms in M". The velocity fields {*{*Qu A 
Z^))") etc are rigid body angular velocity fields in M", where / is position vector 
from the center of the frame, and note that these are also divergence-free. One 
checks that this pairing is also non-degenerate due to the -orthogonality of gra- 
dient fields and divergence-free fields. The Hamiltonian function is again the fluid 
flow kinetic energy written as H{[vl]) := (1/2) J^„{vl + A l^)) A *{vl + 

A /^))) /i, invariant under the change f ^ — )■ t>j + dh. The functional 
derivative is computed as 5H/5[vl] = Vr + A /''))". 

The Lie-Poisson form of Euler's equation in the rotating frame is then given 
by (l50l) with v'" replaced by v^. One obtains 

The equations for the velocity 1-form and the vorticity two-form cUr are then 
given, respectively, by 

+ CM + i*i*^v A I'))) = dg, (51) 

^ + cM^ + d{*i*n,Ai'))) = o, 

In M^, the above vorticity equation reduces to the familiar equation for the vorticity 
vector field of the flow observed in a rotating frame [25 J. Since is invariant in 
time, for the vorticity two-form 

Lu := ujr + d{*{*ily A /'')), 

one readily obtains the evolution equation 

+ A c.) = 0, ^ + £.J, = 0, (52) 

where, as before, cu Aco = /2/i. 

To obtain Ertel's theorem via the above equation in W^, consider a constantly 
rotating frame K3 in M^. This can be viewed as a constantly rotating frame 
in R.^ with the o-axis fixed. The constant angular velocity two-form of K4 is 

Q = Vt^dx Ady — Vtydx Adz + Vt^dy Adz. Consider the vector field f ^ in given 
by 

Vr := {WrJ), 

where Wr is a divergence-free velocity field in and, as before, /(x, y, 2;, t) : 
M'^ —7- M, so that Vr is also divergence-free. One checks that 

Vr + A t')^ := {Wr + A %)% /), 

= {Wr+Xn X /s,/), 

=^ UJ = {(jJr + 2fi, d/), 
^uAuj = {{X^^ + 2Xn) ■ V/) /i, (53) 
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where Z3 is the projection of the position vector / onto M'^, Xq = (fi^, Vly, Viz) and 
the subscript for the parentheses on the right in the first equation denotes that all 
exterior algebraic operations are performed in M'^. 
The system of equations obtained from (|5TI) is then 



dWr f W"^ + P \ 

Euler's equation for Wr in and the conservation law for / are both recovered by 
choosing (7(2;, 2;, o) = (t(;^ + /^)/2 + (Xq x /a) ■ + p/p, where the gradient 
of $ is the centripetal acceleration term in K^. Finally, from (|52|) and (|53l) . one 
obtains Ertel's theorem for the potential vorticity (X^^ + 2Xq) ■ V/ in a rotating 
frame [|25l 

^ ((X^, + 2Xn) ■ V/) = 0, 

where D/DT = d/dt + ■ V. 

A derivation similar in spirit to this is carried out in Subin and Vedan lf30ll 
(section 6) but, seemingly, not with the objective of obtaining Ertel's theorem. The 
authors choose the o-coordinate as time and f = and obtain a final equation 
which is not quite Ertel's theorem. The relation of Ertel's theorem to symmetries 
in Lagrangian and Hamiltonian systems and to symplecticity has been discussed 
previously in the literature — for example, [12611281 1711. See also the recent work of 
Gibbon and Holm |[l2l. 



7 Summary and future directions. 

This paper examines the vortex dynamics of Euler's equations for a constant den- 
sity fluid flow in M^. The focus is on singular Dirac delta distributions of vorticity 
supported on two-dimensional surfaces termed membranes. The self-induced ve- 
locity field of a membrane has a logartihmic divergence. A regularization done via 
the local induction approximation (LIA) then shows that the regularized velocity 
field is proportional to the mean curvature vector field of the membrane rotated 
by 90° in the plane of normals. A Hamiltonian structure for the regularized self- 
induced motion of the membrane is also presented. Finally, the dynamics of the 
four-form a; A a; is examined and it is shown that Ertel's vorticity theorem in M^, 
for the constant density case, is a special case of this dynamics. 

An interesting future direction to pursue is how Hasimoto's transformation 
[lT4l can be generalized and applied to equation (l39l ). and to investigate the re- 
lation of the transformed nonlinear PDF to the nonlinear Schrodinger equation. 
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Constructing membrane models for the dynamics of a; A a; and potential vorticity 
is another direction. As one possibility, models where the potential vorticity is a 
Dirac delta function supported on points could be constructed. In other words, a 
point potential vortex model. However, in general, such models will not provide 
a closed system of evolution equations for the points — unlike for point vortices 
in — since potential vorticity is governed by the evolution of both the vorticity 
field and the passive scalar. Hamiltonian formalisms also exist for incompressible, 
non-constant density and for compressible fluid flows as also for ideal magneto- 
hydrodyanmic flows [[22l[T6ll . Investigating the role of a; A w, or an appropriately 
generalized quantity, in such flows in higher-dimensional spaces is another inter- 
esting direction. 
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A Proof that the strength of a membrane is a well- 
defined constant. 



Proof that Definition (13.11) defines the strength of the membrane. Recall that 
for a filament C in R^, its strength is defined by first considering a de-singularized 
vortex tube. This is a geometric construct and is a tubular neighborhood C C 
of (the curve) C (see, for example, Ana De Silva [|29l for a discussion of tubular 
neighborhoods) such that at every point in dC, igc^ = where igc ■ dC — )■ C 
is the inclusion map. In other words, vortex lines passing through points in dC 
are everywhere tangent to dC. Moreover, it is assumed that there exists another 
tubular neighborhood C D C such that there is no vorticity in C\C. The strength 
of the vortex tube is defined as F := J^i^co, where A C C is again a two- 
dimensional disc-like surface whose boundary dA is a non-contractible loop on 
dC and : ^ — ^ is the inclusion map. Now for any two curves dAi and dA2 
apply Stokes theorem to the part of the tube bounded by Ai, A2 and the part of dC 
that is bounded by dA^. Using the fact that w is a closed form in C and ig^^ = 
delivers the result J^^ i\^uj = J^^ 1^2^- other words, the strength of the tube is 
independent of the choice of dA. This invariance persists as the tube approaches 
the singular limit of the filament. 

For membranes, the procedure is similar. However, some additional steps are 
needed due to the extra spatial dimension involved. By the classical tubular neigh- 
borhood theorem (see again [29]), there exists a tubular neighborhood C C 
of (the surface) S. Analogous to vortex lines being tangent to the boundary of 
a vortex tube in M^, assume that C is such that at each point of dC — a three- 
dimensional 'manifold' — i*g^uj=0, where igc '■ dC — > C is the inclusion. This is 
equivalent to saying that each point of dC there passes a (unique) vortex surface 
S such that i^uj = and S C dC. Again, assume the existence of a neighbor- 
hood C D C such that the vorticity two-form is zero in C\C. Now since C is 
four-dimensional and dco is three-dimensional. Stokes theorem cannot be directly 
applied to C or to any four-dimensional submanifold. 

Consider therefore a tubular neighborhood Bi of a smooth closed curve / G S. 
Bi C C is three-dimensional and dBi C dC. Let is ■ Bi ^ C, then di*^uj = 
since dcu = everywhere in C and the exterior derivative operation on forms 
commutes with puUbacks by mappings i.e di*^uj = i*^duj. Bi is like a vortex tube 
in M'^ but the ambient space is M^. Define the strength of Bi accordingly: 



with A as before (dA a non-contractible loop on dBi) and : A — > is the 
inclusion. As for a vortex tube in R^, the next step is to show that F; is a well- 
defined constant invariant with the choice of OA. Denote by i?i 2 the part of Bi 
with piecewise smooth boundary dBi^2 = AiU A2U dS, where dS is the part of 




(54) 
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dBi bounded by dAi. The next step is to apply Stokes theorem for i^cu in Bi^2- 

f ~ f ~ f ~ 

2 JdS JAi J A2 

where iqs : dS ^ Bi, ia^ : Ai ^ Bi and iA2 : A2 Bi. As shown previously, 
the left hand side is zero. Now, the puUbacks in the integrands on the right are 
the puUbacks by the composition maps o igg etc. which are the same as the 
inclusions igs '■ dS — )• C etc. From the description of dS above, it is obvious that 
igs is the same as igsi ■ dBi — ^ C. And so the three boundary integrals above can 
be written as 



dS JAi 

But, since dBi C dC, igsi is nothing but a restriction of the (previously defined) 
inclusion igc '■ dC — C. It follows that ig^i^ = ^SBj^ac'^ = 0' since i*g(.u = 0, 
where igBi '■ dBi — )• dC is the inclusion. And so the first boundary integral 
vanishes, the second and third represent Ti, as defined by (|54|) . on Ai and A2, 
respectively (keep in mind that the interior of C is everywhere like so that 
and iA2 are really inclusions in M'^) and, after taking into account the orientations 
of the AiS it follows that F/ is the same for dAi and dA2. 

The above shows that the strength of Bi is a constant for a given smooth closed 
curve / G S. The final step is to show that this strength is the same for the Bi 
of any smooth closed curve in E. For this, consider two smooth closed curves 
li,l2 G S which coincide on a finite segment of the curves, denoted by lab, and 
their tubular neighborhoods Bi^ and Bi^, respectively. It follows that Bi^ fl Bi^ 7^ 
0. More importantly, Bi^^ := ^ C Bi^ fl Bi^ is non-empty. Considering a non- 
contractible dA G dBi^^, it is then straightforward to show, using the constancy 
of F^^ and Ti^, that F^^ = Ti^. This geometric construction of two smooth closed 
curves that coincide on a finite segment can obviously also be used to connect any 
two disconnected curves in S — since S is assumed to be arcwise-connected — by 
a third smooth closed curve which shares finite identical segments with each of the 
disconnected curves. A transitivity argument then shows that the Bis of the two 
disconnected curves must have the same strength F^. And this suffices to show 
that the Bi of any smooth closed curve in S has the same strength F/. Passing to 
the singular limit then proves the result ■ 
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